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A selfconsistent calculation of heavy-quark (HQ) and quarkonium properties in the Quark-Gluon 
Plasma (QGP) is conducted to quantify flavor transport and color screening in the medium. The 
main tool is a thermodynamic T-matrix approach to compute HQ and quarkonium spectral functions 
in both scattering and bound-state regimes. The T-matrix, in turn, is employed to calculate HQ 
selfenergies which are implemented into spectral functions beyond the quasiparticle approximation. 
Charmonium spectral functions are used to evaluate eulcidean-time correlation functions which are 
compared to results from thermal lattice QCD. The comparisons are performed in various hadronic 
channels including zero-mode contributions consistently accounting for finite charm-quark width 
effects. The zero modes are closely related to the charm-quark number susceptibility which is also 
compared to existing lattice "data" . Both the susceptibility and the heavy- light quark T-matrix are 
applied to calculate the thermal charm-quark relaxation rate, or, equivalently, the charm diffusion 
constant in the QGP. Implications of our findings in the HQ sector for the viscosity-to-entropy- 
density ratio of the QGP are briefly discussed. 
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I. INTRODUCTION 

In recent years rather remarkable properties of the 
Quark-Gluon Plasma (QGP) have been discovered, no- 
tably a liquid-like behavior deduced from the collective 
expansion of the medium formed in collisions of heavy 
nuclei at the Relativistic Heavy Ion Collider (RHIC). Hy- 
drodynamic expansion models, based on the assumption 
of local thermal equilibrium, require very short equili- 
bration times, of the order of Ttherm — Ifm/c or less, 
to account for the experimental findings. The micro- 
scopic origin of the rapid thcrmalization remains one of 
the outstanding problems raised by the RHIC data. It 
is further complicated by the fact that the very notion 
of equilibrium erases information about its origin. A 
phenomenologically suitable probe of thcrmalization in 
heavy-ion collisions should therefore interact strongly at 
the thermal scale but with a relaxation time on the or- 
der of the system's lifetime (at RHIC tqgp — 5fm/c). 
The logical choice are thus heavy quarks, i.e., charm and 
bottom {Q = c,b), for which thcrmalization is expected 
to be delayed by a factor mq/T, quite comparable to 
TQGp/nherm (scc, e.g., Rcf. [l| for a recent review). 

Bound states of heavy quarks (charmonia and bot- 
tomonia) are believed to probe the QGP from a some- 
what different angle [H-Q. The dissolution pattern of 
the quarkonium spectrum in matter is possibly the most 
direct way of diagnosing color-Debye screening of the 
basic QCD force, fgg = ^^^qq('')j as a function of 
temperature (and/or density). In-medium quarkonium 
spectroscopy therefore reveals insights into the dcconfin- 
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ing hadron-to-parton phase transition in QCD. In par- 
ticular, the use of potential models has been revived 
recently, largely triggered by the prospect that an in- 
medium potential can be extracted model-independently 
from lattice QCD (IQCD) computations. Furthermore, 
IQCD generates euclidean-time (r) correlation functions 
of quarkonia with good accuracy in different hadronic 
channels, which provide useful constraints on model cal- 
culations of spectral functions in the timelike regime of 
physical excitations. The extraction of reliable informa- 
tion from such comparisons requires a realistic modeling 
of the continuum part (scattering regime) of the quarko- 
nium spectral functions, i.e., not only its bound-state 
part. The account of interactions in the near-threshold 
region is particularly important to describe situations 
were bound states dissolve into the continuum, as in 
the case at hand. In connection with correlator analy- 
ses, a comprehensive treatment of bound-state and con- 
tinuum regimes has been performed using Schrodinger 
phase shifts 0] , a thermodynamic T-matrix Q or a non- 
relativistic Green's function Q approach. Another im- 
portant ingredient to a realistic description of quarko- 
nium spectral functions in the QGP are medium ef- 
fects in the single-particle properties, i.e., in the heavy- 
quark (HQ) propagation. These are closely related to 
HQ transport [10, and the so-called zero-mode con- 
tributions to quarkonium spectral and correlation func- 
tions [H], [l^. In particular, finite- width effects on both 
quarkonia and heavy quarks have received little attention 
thus far 

In the present paper we further develop our previous 
study of heavy quarks and quarkonia in the QGP using a 
thermodynamic T-matrix approach [loj . First, we go be- 
yond the quasi-particle approximation for heavy quarks 
and perform a selfconsistent calculation of the HQ self- 
energy and the in-mcdium heavy-light T-Matrix. The 
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improved HQ spectral functions are then implemented 
into an off-shell calculation of the quarkonium spectral 
functions. Second, we expand our comparison with Eu- 
clidean correlator ratios to the scalar, vector and axi- 
alvector channels. In these channels the presence of zero 
modes [11113, resulting from scattering of a single heavy 
quark on medium constituents ( "particle- hole" excita- 
tions), is known to impact the large-r behavior of the 
in-medium correlators appreciably. Thus far, these con- 
tributions have been estimated in quasiparticle approxi- 
mation neglecting finite width effects [3, [H, [H, [IE [111 ■ 
In a treatment of heavy quarks consistent with the T- 
Matrix we evaluate the zero-mode contribution account- 
ing for their full spectral function, in particular finite- 
width effects. This automatically yields predictions for 
the charm-quark number susceptibility, which can be 
compared to lattice data and enable a more reliable 
extraction of in-medium quasiparticle masses for heavy 
quarks. 

The paper is organized as follows: In Sec. [Tl] we give 
a short overview of the T-Matrix formalism employed 
in our calculations, specifically discussing its extension 
to off-shell dynamics in the charm-quark propagation. 
Section IIIII is devoted to properties of cahrmonia in the 
QGP, i.e., their spectral functions including zero modes 
fSec. lIIl"A^ . cuclidean correlator ratios in various mesonic 
channels (Sec. IIII B[) and a comparison of numerical re- 
sults using two different input potentials to IQCD data 
(Sec. IIII C|) . In Sec. IIVI we first analyze the selfconsis- 
tently calculated charm-quark selfenergies (Sec. IIV A[) 
and then apply those to obtain the thermal relaxation 
rate and a schematic estimate of r]/s (Sec. IIVB|) . as well 
as the charm-quark number susceptibility (Sec. IIV Cp . 
We conclude in Sec. |Vl 



II. OFF-SHELL T-MATRIX AT FINITE 
TEMPERATURE 

The formalism for calculating the T-Matrix for quark- 
quark and quark- antiquark scattering and bound states 
in the QGP, using two-body potentials estimated from 
heavy-quark free energies computed in lattice QCD 
(IQCD), has been developed in Refs. [i,[IS[Il. Here we 
recollect the main elements while referring to Ref. [I3| for 
further details (e.g., a discussion of relativistic corrections 
and constraints from vacuum spectroscopy and the high- 
energy perturbative limit). Starting point is the Bethe- 
Salpeter Equation which after a 3-dimensional (3D) re- 
duction and partial- wave expansion turns into a ID inte- 
gral equation for the scattering amplitude (T-Matrix), 

00 

Ti^a{E;q\q)^Vi,a{q',q)+'^ j dk ViAq' , k) 



xGi2{E;k)Ti,a{E;k,q) , 

(1) 



in a given color channel (a) and partial wave (l); the rel- 
ative 3-momentum moduli of the initial, final and inter- 
mediate 2-particle state are denoted by g = |q |, = |q' | 
and k = |k |, respectively (we restrict ourselves to vanish- 
ing total 3-momentum, P = 0, of the 2-particle system). 
The explicit form of the intermediate 2-particle propa- 
gator, G12, depends on the 3D reduction scheme |18l - 
[2]| (in slight deviation from Refs. [^ [13, [13 we here 
absorb the Pauli blocking factor, (1 — 2fp), into G12). 
In the following, we focus on the Thompson scheme, 
since the Blankcnbeclcr-Sugar scheme was found to gen- 
erate some overbinding in the vacuum quarkonium spec- 
trum [lo| . The T-matrix will be applied in both heavy- 
heavy and heavy-light quark channels for which a static 
(potential) approximation can be justified (i.e., the en- 
ergy transfer is paranicterically suppressed compared to 
the 3-momentum transfer, Aqo {Aqj'^/mq << Aq, 
where Aq = q' — q). For the two-body potential, 
Vz,a, we follow our previous work [l^ using either the 
heavy-quark free (F) or internal (U) energy computed 
in IQCD, implemented into a field-theoretical model for 
color-Coulomb and confining terms |22| with relativistic 
corrections (e.g., the Breit interaction |23l - |25| to account 
for color-magnetic effects). To ensure the convergence 
of the Fourier transform from coordinate to momentum 
space we subtract the infinite-distance limit of the po- 
tential according to 

Va{r;T) = Xair,T)-X{r ^ ^,T) , X = F or U , (2) 

and interpret XociT)/2 = X{r -s- cx),T)/2 as a 
temperature-dependent correction to the bare HQ mass 
(real part of selfenergy) induced by the condensate asso- 
ciated with the confining force term. The use of cither U 
or F as potential is believed to bracket the uncertainties 
in this identification. 



A. 2-Particle Propagator 

In all previous applications of potential models in the 
QGP the quarks have been treated as quasiparticles with 
either vanishing or constant [1, [I^ width. However, as 
pointed out in Ref. [lol |. this approximation implies that 
quarkonium spectral functions, Pa{E) [a: quantum num- 
bers of the composite mesonic (or diquark) state) , do not 
possess the proper low-energy limit, pa{E 0) (x E. 
While this is not expected to significantly impact the 
mass and binding energy of the bound states (for which 
the total energy is of order twice the HQ mass), the cu- 
clidean correlators calculated below do actually involve 
an integration over Pa{E) starting from E = Q with a 
thermal weight which diverges for — > 0. In previ- 
ous studies this problem has been evaded by introduc- 
ing a low-energy cutoff on the spectral function (suf- 
ficiently below the lowest bound state as to not affect 
the correlator). However, recent studies of quark spec- 
tral functions (26| find considerable low-energy strength 
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in quarkonium spectral functions, e.g., due to particle- 
hole like structures in the HQ propagator. Thus a more 
elaborate treatment in the T-Matrix equation, properly 
accounting for off-shell dynamics, is in order. 

The general off-shell expression of the uncorrelatcd 2- 
fermion propagator at finite temperature, 6*12, figuring 
into the T- matrix equation, can be derived, e.g., within 
the Matsubara formalism. One has 

Gi2{nx,k) =TY,Gi{z,,k)G2{nx- z,,-k) (3) 
where 



G.(w,k) = 



cj^ — fc^ — m~ — 2miT,i{uj, k) 



(4) 



(i = 1, 2) denotes the (scalar part of) the 1-particle prop- 
agator, and Zn = i{2v+\)'KT (SIa) arc fcrmionic (bosonic) 
Matsubara frequencies. Using the spectral representa- 
tions of the in-medium (retarded) single-particle propa- 
gators, 



doj' Pi{i-u' , k) 



= -2 ImG. , (5) 



the Matsubara summation in Eq. ([3]) can be performed 
explicitly. After analytic continuation to the real axis 
(r^A E), the positive-energy contributions to the 2- 
particle propagator take the form 



X mi 1712 



E — uj — uj' + ie 



(6) 



to gT with a coefficient to approximately match IQCD 
computations of the energy density of the QGP in quasi- 
particle approximation. In addition, we supplement a 
moderate imaginary part for the light-quark selfenergy, 
Im Eg = —0.05 GeV, which is in the range of values ob- 
tained in Ref. [iTj in a similar approach for light quarks, 
and comparable to what we obtain for heavy quarks (our 
results for quarkonium spectral function and HQ trans- 
port are largely insensitive to the light-quark width) . 

The evaluation of the Matsubara sum in Eq. ([3|) also 
generates a contribution that solely arises from thermal 
excitations (see, e.g., Ref. H^]). It involves one positive- 
and one negative-energy part of the quark spectral func- 
tions which is why we referred to it as a "particle-hole" 
contribution above. For mcsonic channels of the T- 
matrix, 12 — QQ, its imaginary part can be cast into 
the form 

Im G{^{E, k)^-J ^Pq(^, k)p+{E + lo, k) 

X [/«(L.,fc)-/Q(ii; + c^,fc)] . (9) 

Here, the negative-energy part of the Q spectral func- 
tion has been turned into the positive-energy part of the 
Q spectral function, i.e., the Q line in the original QQ 
propagator has been "turned around" . Together with 
the Fermi distributions, the interpretation of this contri- 
bution to 6*12 becomes apparent: an incoming Qj pre- 
existing in the heat bath according to /*^(w, fc), scatters 
into a Q with energy E + lo, Pauli-blockcd according to 
f'^{E + oj,k). G12 is, in fact, precisely the zero-mode 
contribution discussed in the context of quarkonium cor- 
relators m, [13, [iH. We return to its evaluation in 
Sec. HEl below. 



where 



-Im 



(fc) w - Wi(fc) - I]i(w, fc) 



,(fc) = Jk^ + 



(7) 



denotes the positive-energy part of the quark spectral 
functions. The factor mim2 in Eq. ([6]) is specific to the 
Thompson reduction scheme, rendering the appropriate 
expression for G12 in the limit of on-shell quarks [2]| . 
For the in-medium HQ mass, as mentioned above (recall 
Eq. the infinite-distance value of the HQ potential 
is identified with a "mean-field" contribution, i.e., as a 
real part of a selfenergy. 



'^ + E|^Mr) , S^^^(r)^X^(T)/2. (8) 



The additional selfenergy contribution figuring into the 
quark propagator in Eq. (O is generated from interac- 
tions with thermal light quarks and antiquarks and is 
therefore distinct from Xoo{T). It will be computed 
from the heavy-light quark T-matrix, as discussed below. 
For calculating the latter, we need to specify light-quark 
masses which we model by a thermal mass proportional 



B. Single-Quark Selfenergy 

The HQ selfenergy, Y,q(uj, k), due to interactions with 
light anti- / quarks in the heat bath is calculated from the 
heavy- light (Qq) T-Matrix by closing its light-quark line 
with a quark propagator weighted with a Fermi distribu- 
tion function. Using the Matsubara formalism one can 
express its imaginary part as [TtI . [27j 



-Emin 

X ImMqQ {E,LL>,LL>q{p),k,p) 
x[fp{u,ip))+fBiLO+UJ,{p))]{10) 



with 



MgQ{E,uj,uj\k,p) 



T-qrriQ 



Wg((7cm) WQ((/cm) 

47r da[TaAE,qcm) + 5TaAE,qcm)] (11) 



a— 1,1 
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where 



9cm(^.fc(4), 
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The spin-isospin factor dsi = 4:Nf counts the degener- 
acy of available meson (or diquark) states, e.g., with total 
spin-0 and -1 for 5'- wave heavy-light scattering [l3l. In 
the expression for the HQ selfcncrgy, Eq. ([TU]). the ther- 
mal light quarks are treated as zero- width quasi particles 
so that their spectral functions can be replaced by S func- 
tions (but with thermal mass ^ gTy. We recall that the 
gluonic contributions are entirely attributed to the mean- 
field type condensate term (for HQ transport, we include 
heavy-light interactions to leading order in perturbation 
theory which does not generate an imaginary part in the 
scattering amplitude). The real part of the selfenergy 
is obtained by a dispersion relation which is a preferred 
procedure in a selfconsistent prescription since the nor- 
malization of the spectral functions can easily be guar- 
anteed. Our framework is similar to the one utilized in 
Ref. [l3| in the light-quark sector. However, in there the 
calculations of the T-matrix were restricted to on-shell 
selfenergies while here we use the full off-shell selfenergy 
of the heavy quark which, in particular, enables to es- 
tablish the correct low-energy behavior of the mcsonic 
spectral functions (in addition to refinements in the im- 
plementation of the potential as developed in Ref. [l3| the 
non-potential corrections are expected to be significantly 
larger for light-quark interactions). 

Let us finally comment on the issue of imaginary parts 
in the potential and selfenergy. Using effective field the- 
ories at finite temperature it has been found that the 
two-body potential operator figuring into a Schrodinger 
equation acquires an imaginary part [28l - l3l| . Diagram- 
matically, this implies that the potential possesses on- 
shcU cuts corresponding to dissociation processes of the 
composite (bound) state. Let us inspect two commonly 
discussed ( "leading-order" ) processes for the case of char- 
monium bound states. For large binding energy, Eb > T, 
the dominant process is gluon dissociation, g + Jl ip ~^ 
c + c, first analyzed by Bhanot and Peskin [32|. In 
the language of effective field theory, this corresponds 
to the so-called color-singlet to -octet transition mech- 
anism (up to final-state interactions in the octet chan- 
nel, which, however, are suppressed by l/N^ and thus 
numerically negligible). In our T-matrix formalism, the 
inclusion of this process would require a coupled-channel 
treatment, with a cc-gluon intermediate state, whose cut 



^ This is mainly done for numerical reasons. We have checked that 
using an off-shell quark spectral function leads to very similar 
results. 



exactly produces the corresponding decay channel. Such 
a calculation is beyond the scope of the present work. 
For small binding energies, Eb < T (as relevant for ex- 
cited states or the case of reduced in-medium J/ip bind- 
ing energies close to its dissolution temperature), gluo- 
dissociation is phase-space suppressed and thus super- 
seded by the Landau-damping phenomenon in the (space- 
like) onc-gluon exchange of the potential. This is pre- 
cisely the imaginary part of the two-body potential dis- 
cussed in Refs. [H, [s^, HH. Physically, it corresponds 
to "quasifree" dissociation [s^, p + J/tp p + c + c, 
induced by thermal partons p = g,q,q (which generate 
the imaginary part of the one-loop correction to the ex- 
changed gluon). In the limit of small binding (or large 
charmonium size), the incoming thermal partons with 
energy '^T do not sense the size of the bound state and 
thus the scattering effectively happens on the level of an 
individual charm (or anti-charm) quark. Thus, in the 
T-matrix formalism this process is encoded in the selfen- 
ergy of a single (anti-) charm quark, which is included in 
our calculations. For large binding, quasifree dissociation 
ceases since the thermal parton does not resolve the sub- 
structure of the color neutral bound state [1^ [s^ . On 
the other hand, when formulated as a potential contribu- 
tion, the large-distance limit of its imaginary part coin- 
cides with twice the imaginary part of the charm-quark 
selfenergy [2^. Therefore our evaluation of the correla- 
tors within the T-Matrix approach using a real potential 
is well in line with the use of a complex potential in a 
Schrodinger treatment. In the T-matrix, imaginary parts 
are generated via unitarization of intermediate (on-shell) 
states including single-particle selfenergy contributions. 



III. QUARKONIA 

In this section we apply the selfconsistent T-matrix to 
charmonia with special consideration of the zero-mode 
contribution (Sec. IIII A[) to the Euclidean correlators 
in different mcsonic channels, a = S,PS,V,AV, i.e., 
scalar, pscudoscalar, vector and axialvector, respectively 
fSec. IIII B|) . followed by a discussion of numerical results 
using either C/ or T as underlying two-body potential 

(Sec. mini). 



A. Spectral Functions and Zero Modes 

With the cc T-Matrix from Eq. ^ we proceed to de- 
termine the charmonium spectral function including both 
bound and scattering states as in Ref. [l3l . The T-matrix 
signifies the rescattcring contribution to the correlation 
function which is schematically given by Q 



G = G° + G°TG° 



(13) 
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where G° denotes the free 2-particle loop. This can be 
be represented diagrammatically as 



G = 




(14) 



where compared to Eq. © the loop now also includes an 
integration over relative momentum, k, as well as ver- 
tices (denoted by dots) specifying the quantum-number 
channel a. The spectral function is then defined as usual 
by p = -2ImG. 

Quantitative comparisons to euclidean correlator ra- 
tios as "measured" in IQCD lSil. 13511 require the inclusion 
of zero-mode contributions [121 Il3| which turn out to be 
different for different meson channels a. This, in partic- 
ular, lifts the spin degeneracy within 5-wavc {PS — V) 
and i-*-wave (S — AV) states of the T-matrix. This is not 
surprising since HQ symmetry is not expected to be valid 
for zero modes. A rclativistic evaluation of the vertices 
figuring into p™ = — 2Im is therefore in order. To be 
consistent with the treatment of the T-Matrix we evalu- 
ate this contribution beyond the zero- width quasiparticle 
approximation (as applied in the literature to date) by 
taking into account the finite width of the quark. As in 
previous literature, we focus on the leading part of the 
zero-mode contribution which does not involve a two- 
body potential contribution^. Augmenting Eq. (|9]) with 



relativistic vertices one obtains jl2| 



pT{P) = 2N,j ^Tr[(^ + me)r(/ + m,)rt] 
xpt{k)pt{r) [./^(fc°)-/^(r°)] 



2N, 



i2ny 



■4:Fr{k,r, rric 



xpt{k)pt{r) [/^(fc°)-/^(r°)] (15) 
where we introduced the notation 



r = P + k 



r e {1,75 ,7o,7»,757j (16) 



with r, k and P = {E, P) denoting 4- vectors for the 
purposes of the above two and the following equation 
only. The kinematic factors arising from the different 
Dirac structures take the form 

Fi(fc,r, m) = k ■ r + iTL^ 
Fj,^{k,r,m) = k ■ r — m\ 
Fyg{k,r,m) = fco • ro + k • r -|- 



Fj. (fc, r, m) 
^j^^.{k,r,m) 



3r°A:°-r-k 
3r°/c° - r -k- 



3m^ 



(17) 



a summation over spatial indices i is implied for the vec- 
tor and axialvector case. In the previously adopted zero- 
width approximation the zero-mode spectral function can 
be simplified to [l^ 



pr(^,P^0) = 2^i?5(i?)xc(T) 



(18) 



where 



2Nr 



(27r)3 



Cl + C2 - 



5/^K(k)) 
Wc(k)2y dujc{^) 



(19) 

denotes a generalized susceptibility with coefficients cor- 
responding to different quantum number channels ac- 
cording to 

r = l => ci,2 = 2,-2 r = 75 =^ ci,2 = 0,0 

r = 7o ^ci.2 = 2,0 r = 7, ^ ci,2 =0,2 

r = 757i ^ ci.2 =6,-4 . 

In our numerical calculations reported below we evalu- 
ate p™ directly from Eq. ([T5|) using the positive-energy 
charm-quark spectral function, Eq. ([7|), with the selfcon- 
sistcntly determined off-shell sclfencrgy. This puts the 
treatment of the zero-mode contribution on the same 
level as the cc scattering and bound-state part using the 
2-particle propagator G12 in Eq. ©. 



B. Euclidean Correlator Ratios 

The transformation of the spectral function to the eu- 
clidean correlator is given by 



G„(r,T) 
ICir,E,T) 



dE 
2^ 

cosh [E{t ~ l/2r)] 



p^{E,T)+p^^{E,T)] IC{t,E,T) 

(20) 



sinh[£;/2T] 



As before, we focus an vanishing total 3-momentum, 
P = 0, of the QQ. The temperature kernel /C imprints 
an exponential r dependence on the correlator. To bet- 
ter exhibit the in- medium modifications of Gq (t, T) it is 
therefore common to analyze the correlator ratio 



i?g(r,T) 



G„(T,r) 

G-(r,T) 



(21) 



where the denominator is the so-called reconstructed cor- 
relator which is evaluated with the kernel JC at tempera- 
ture T but with a "reference" spectral function typically 
taken as the vacuum one (or at low temperature T* where 
no significant medium modifications arc expected). 



^ A non-perturbative treatment of the particle-hole interac- 
tion would require a relativistic treatment solving the Bethe- 
Salpeter-Equation which goes beyond the potential approxima- 
tion adopted here. 



Numerical Results 



Before turning to the numerical results for the finitc-T 
spectral functions and correlator ratios let us briefly sum- 
marize our input quantities, largely as given in Ref. (lo| 
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FIG. 1: (Color online) Charmonium spectral functions in the pseudoscalar (upper left panel) and scalar (upper right) channel 
for different temperatures using the internal energy, [/, as QQ potential. In the middle and lower rows we show the corre- 
sponding euclidean correlator ratios in the pseudoscalar (middle left), scalar (middle right), vector (lower left) and axialvector 
channels (lower right), where the latter three include zero-mode contribution. The comparison to the lattice data (indicated 
by symbols) [s^ is made in absolute units of r (limited to XjlT at each temperature), with Tc — 210 MeV for the lattice data 
and Tc ~ 196 MeV underlying the potential for the T-matrix calculations (e.g., 1.4 T^^^ ~ 1.5 T^^™'''). 



(where more details can be found, including extensive 
analysis of the associated uncertainties). Our poten- 
tials (F or [/) are based on fits to the lattice results for 
(2+l)-flavor QCD from Refs. [1113 ("potential 1" in 



Ref. [T3l)- In vacuum our ofF-shell calculations reproduce 
the results of Ref. [lo| since in the limit of vanishing c- 
quark width the propagator 6*12 reduces to the standard 
Thompson form. The bare charm-quark mass is set to 
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FIG. 2: (Color online) Same as Fig. [T] but using the free energy, F, as QQ potential. 



771° = 1.264 GeV, which, together with a vacuum selfen- 
ergy of = 0.6 GeV, gives a fair description of the 

spin-average masses of J/ip-r]c, ^' and Xc states. 

The full off-shell treatment with more realistic c-quark 
propagators in the present work leads to moderate but 
significant changes for the in-medium results. In ad- 
dition, the zero-mode contributions in the S and 
AV channels (recall Eq. p^ *) have a marked impact on 
the corresponding correlator ratios. For a comparison 



to IQCD correlator ratios we choose the results of the 
Nf = 2 computations from Ref. [sH, which correspond 
more closely to our input potentials than quenched cal- 
culations. The critical temperature in the simulations of 
the IQCD correlator ratios [H is ~ 210 MeV, which 
is not far from the one underlying our potential [36j . 
Tc ~ 196 MeV. Therefore, rather than normalizing to 
the different Tc's, we compare our results for the cor- 
relator ratios with the IQCD data in absolute units of 
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euclidcan time, r. Each correlator ratio is plotted up to 
the midpoint, r = 1/2T, of the total r range (the correla- 
tors are symmetric about this point). Thus, comparable 
temperatures are easily identified by the same r range in 
the plotted ratio (if one were to normalize all results to 
Tc, IQCD data at given absolute temperature would be 
compared to T-matrix calculations at slightly lower T). 

The in- medium results using the internal energy, U , 
as potential are summarized in Fig. [l]in terms of the S*- 
and P-wave spectral functions in the upper panels (de- 
generate for rjc-J/tlj and for Xc states) and the correlator 
ratios in the middle and lower panels (where the degen- 
eracies arc lifted by the zero modes). Since the (magni- 
tude of the) imaginary part of the selfconsistently calcu- 
lated c-quark selfenergy turns out to be around 0.050- 
0.100 GeV for low-momentum on-shell charm quarks 
(cf . Fig. [3] below) , the closest comparison to our earlier 
quasiparticle results is for the case of a constant (en- 
ergy and 3-momentum independent) imaginary part of 
ImSg = -0.05 GeV (Fig. 14 in Rcf. Hfl]). In the S- 
wave spectral function we find slightly more attraction 
for the ground-state peak for temperatures below 1.5 Tc 
(for higher T it dissolves, as in Ref. [i3). Despite the 
larger on-shell width in our present treatment, which is 
up to twice as large as in our previous quasiparticle cal- 
culations, the width of the J/'^-rjc peaks is very simi- 
lar. This is a direct consequence of the energy and mo- 
mentum dependence of the selfenergy which decreases 
considerably off-shell (cf. Fig. [3] below) and thus reduces 
the "operative" width of c-quarks in cc bound states. It 
becomes apparent in the rjc correlator ratio (where no 
zero-mode is active) which at the lowest considered tem- 
perature (1.2 Tc) drops to about 0.85 compared to 0.9 in 
Ref. even though the peak position of the bound 
state is shifted to slightly lower energies in the present 
calculation (which tends to increase the large-r correlator 
ratio). This finding shows that a proper off-shell treat- 
ment is warranted to correctly account for the low-energy 
strength in the spectral function which has significant im- 
pact on the correlator ratio at large r. We are not too 
concerned that the drop to 0.85 is noticeably larger than 
in the IQCD data since we have neglected several effects 
which will contribute further low-energy strength to the 
spectral function, e.g., imaginary parts in the c-quark 
selfenergy from scattering off thermal gluons, or coupled 
channels in the cc T-matrix such as DD and ccg (in- 
ducing singlet-to-octet transitions). These are expected 
to be especially relevant close to Tc (Z?-mcson states will 
form close to Tc and gluo-dissociation is efficient for large 
charmonium binding, Eb > T). On the other hand, we 
note that the PS correlator ratio is remarkably indepen- 
dent of temperature and closer to one for intermediate 
and small r than in Ref. [l^l (e.g., no more than ^2% 
above one), which improves the agreement with IQCD. 
In the V, S and AV channels the zero-mode contribu- 
tions lead to a marked enhancement of the correlator ra- 
tios at large r, especially for the P-wave channels where 
the ci coefficient is non-zero (for the vector channel, we 



sum over the spatial components only, corresponding to 
r = 7i in Eq. ([15]). Compared to a zero-width treat- 
ment of the zero-mode contribution, corresponding to 
Eq. (jl9p . the inclusion of a finite quark width increases 
the correlator ratios at large r by ca. 0.1-0.15. Over- 
all, the agreement of the the calculated correlator ratios 
with Nf = 2 IQCD data (sEj included in the middle and 
lower panels of Fig. [1] is fair. The largest discrepancies of 
ca. -^30% at both intermediate and large r occur in the 
AV (xci) channel where the zero- mode contribution is 
the strongest. The extracted melting temperatures of the 
charmonium bound states are close to our earlier determi- 
nation with finite but constant width [lo| . i.e., ca. 1.5 Tc 
for the J/* and below 1.2 Tc for ah other states (\1/', Xc)- 
The results obtained using the free energy, P, as poten- 
tial are compiled in Fig. [5] As in previous work j9l. [l0l . l39| . 
we find a much stronger suppression of the bound states 
compared to using [/, with a melting temperature for 
the S'-wave ground state of 1.2 Tc or even lower, whereas 
the P-wave spectral function is already structureless at 
this temperature. Nevertheless, the correlator ratio in 
the ?7c channel (no zero mode) is surprisingly T-stable 
and close to one: the loss of low-energy strength due to 
the dissolved bound state is compensated by a reduced 
cc threshold in connection with a nonperturbative rescat- 
tering strength in the threshold region generated through 
the T-matrix Q. Indeed, the in- medium charm-quark 
mass (correction) following from Eq. ((8)) for Poo is signif- 
icantly smaller than for Uao- This, however, generates a 
markedly enhanced zero-mode contribution to the corre- 
lator ratios in the V , S and AV channels, compared to 
the case with U as potential (the effect of the selfcon- 
sistently calculated finite quark width in the zero modes 
is a ca. 10% enhancement at large t). In particular in 
the P-wave channels, the large-r is quite a bit off the 
Nf~2 IQCD data. Apparently a potential closer to the 
internal energy is favored (with larger m* but stronger 
cc binding). 

IV. OPEN-CHARM TRANSPORT 

We now turn to examining the numerical results for 
the single charm-quark properties in the QGP, specifi- 
cally its selfenergy (mass correction and scattering rate; 
Sec. lIV A| . thermal relaxation rate fSec. lIVB]) and num- 
ber susceptibility (Sec. IIV C|) . We recall that these are 
selfconsistently computed via numerical iteration with 
the T-matrix in all heavy-light quark channels, Eq. ([l]). 

A. Selfenergy 

The real and imaginary parts of the charm-quark self- 
energy following from the T-Matrix are displayed in 
Fig. 13] as a function of quark energy for vanishing 3- 
momcntum and for three temperatures. The general 
structure is that of a maximum around the on-shell en- 
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FIG. 3: (Color online) Real and imaginary parts of the charm-quark selfenergy at vanishing 3-momentum in a QGP at different 
temperatures for the scenarios where the internal (left panel) or free energy (right panel) is employed as potential. The 
imaginary parts are negative definite while the real parts change sign close to the vertical lines which indicate the effective 
charm-quark mass (on-shell energy) at the respective temperature. 



ergy, uj ~ m*, in the imaginary part, associated with a 
typical zero crossing in the real part. The peak structure 
is caused by Feshbach type resonances (or threshold en- 
hancements) around the heavy-light quark threshold in 
the T- matrix These features are the reason that, 

upon (off-shell) integration over w in the two-particle 
propagator (T- matrix), the contributions from the real 
part largely cancel while the width effects drop signifi- 
cantly faster (and thus on average are smaller) than for a 
fixed (constant) quasiparticle particle width as employed 
in previous work [1, [13, [13 • 

More quantitatively, in the case of U as potential (left 
panel in Fig. |3|), the on-shell width of the charm quark 
reaches a value of up to Fc = —2lmSc — 250 MeV at 
the lowest considered temperature (1.2 JV), quite simi- 
lar to what has been calculated in Ref. [IJ] using differ- 
ent IQCD inputs for U. The magnitude decreases with 
increasing temperature (not as much as in Ref. [l3|). 
which is quite remarkable given the fact that the ther- 
mal densities of the thermal anti-/quarks decrease ap- 
preciably (in the "T-p" approximation, the selfenergy 
would be directly proportional to the light-quark den- 
sity, Ec ^ TcqUq). It clearly reflects the weakening of 
the two-body interaction as temperature increases, due 
to color-charge screening in both the Coulomb and con- 
fining parts of the potential. In other words, the maximal 
interaction strength is realized at the lowest temperature, 
just above Tc. 

For the .F-potential the c-quark selfenergy is reduced 
by about a factor of 6 at the lowest temperature, which 
reduces to a factor of '-^2-3 at 2Tc. In this scenario, no 
Feshbach resonances form, but rescattering embodied in 
the T-matrix still produces a notable threshold enhance- 
ment which induces an on-shell width of ca. 40-70 McV. 



B. Thermal Relsixation Rate 

Next we turn to the calculation of the thermal relax- 
ation rate of charm quarks in the QGP, employing a 
Fokker-Planck treatment following Ref. [40]. The rela- 
tion of the drag coefficient, A{p), to the T-matrix has 
been elaborated in Ref. [iJli including all color configu- 
rations in c-q and c-q scattering in S- and P-waves. For 
a more realistic evaluation of the total coefficient, we add 
to the T-matrix contribution the effect of c-quark scat- 
tering off thermal giuons which we approximate with the 
leading order perturbative diagrams (using as=0.4) in- 
cluding Debye screening masses in the exchange propaga- 
tors (t-channel gluon exchange gives the dominant con- 
tribution). In Fig. U the full results for A{p) are com- 
pared to perturbative calculations in which scattering off 
both anti-/quarks and giuons is obtained from the LO 
diagrams. Compared to the quasiparticle treatment in 
Ref. [l^, the full off-shell treatment leads to a ca. 10% 
increase of the drag coefficient at low 3-momenta, while 
the deviations are small at momenta of fc > 2 GeV. This 
is well in line with the finding in Ref. [lo| that the depen- 
dence of the drag coefficient on the 3D reduction scheme 
of the T-matrix equation (relating to its off energy-shell 
behavior via different 2-particle propagators) is small. 
Thus, we confirm as a robust feature that c-quark ther- 
malization using the nonperturbative T-matrix is accel- 
erated over pQCD calculations by about a factor of 4 (2) 
when using U (F) as potential. 

In addition to the direct calculation of the relax- 
ation rate following from the collision term in the Boltz- 
mann (or Fokker-Planck) equation with the heavy-light 
T-matrix [l^, [T^] , the off-shell framework set up in the 
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FIG. 4: (Color online) Charm-quark relaxation rate as a function of 3-momentum for different temperatures using either the 
T-Matrix plus pQCD gluon scattering [lol ] (upper 3 lines in left panel) or a pQCD calculation for scattering off gluons and 
anti/quarks (lower 3 lines in left panel). The T-matrix is computed using U (left) or F (right) as potential while pQCD 
contributions are obtained with Qs=0.4. 
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TABLE I: Comparison of charm-quark relaxation rates, 
7c(fc = 0;T), as obtained from the heavy-light T-Matrix [l^ 
(2. column) and from a Kubo formula for the diffusion coef- 
ficient (3. column) using the low-energy limit of the heavy- 
quarkonium spectral function, Eq. (|22p . for the two different 
scenarios for potential (upper and lower table). The last col- 
umn lists a schematic estimate of the viscosity to entropy- 
density ratio using the kinetic-theory relation, Eq. p3p . in- 
cluding the contributions from scattering off thermal gluons 
(the quoted ranges only reflect the variation in the 2 preceding 
columns; the true uncertainty is larger). 



present paper enables an alternative method, namely 
from the zero-energy limit of the heavy-quarkonium spec- 
tral function [ill EH j 

7. = ^ i^.^^hm^^ (22) 



(Xc = Xoo: charm-quark number susceptibility, Dg'- spa- 
tial diffusion constant). The values extracted from this 
method correspond to the zero-momentum limit and are 
compared to the pertinent values from the collision in- 
tegral using the heavy-light T-matrix in Tab. |l] (for a 
consistent comparison the contribution from scattering 
off thermal gluons is not included). Note that the latter 
involves the Qq scattering amplitude squared, while the 
imaginary part of the quarkonium spectral function is 
basically determined by the imaginary part of the single- 
quark spectral function which, in turn, is obtained from 
the imaginary part of the heavy-light scattering ampli- 
tude, cf. Eq. (|10|) . Wc find that the extraction from the 
low-energy limit of the charmonium spectral function in 
the vector channel tends to give larger (smaller ) values 
at lower (higher) for both [/ or F as two-body poten- 
tial. One of the uncertainties which is presumably re- 
flected by these deviations is the Fokkcr-Planck approx- 
imation when using the heavy-light T-matrix to evalu- 
ate the diffusion coefficient. For example, in Ref. [i^l 
it has been found that a T)-meson resonance model in 
the QGP can lead to a violation of the Einstein rela- 
tion by underestimating of 7c by up to 10% at temper- 
atures of T ~ 300 MeV, while the agreement is closer at 
lower temperatures (less than 5% for T < 250 MeV; note 
that the T'-meson like correlations in the T-matrix dis- 
solve for temperatures above 1.5 Tc ~ 300 MeV). Overall, 
the discrepancies between the two methods are within 
^^30%, which, given the different schemes of evaluating 
this transport coefficient, is not too bad. The systematic 
trends are similar when using F or U as potential, and 
the large difference between these 2 scenarios is robust. 

Let us end this section by a schematic evaluation of 
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FIG. 5; (Color online) Charm-quark number susceptibility as a function of temperature normalized to the zero mass limit, 
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the widely discussed ratio of viscosity to entropy density, 
77/s. Using kinetic theory, one can roughly relate this 
quantity to the spatial HQ diffusion coefBcient as 0] 



s 5 



(23) 



The numerical coefficient probably constitutes a lower 
limit, applicable to a weakly coupled system; it is most 
likely larger for strongly coupled liquids, e.g. '^1/2 in 
AdS/CFT. Nevertheless, the results for U suggest that 
77/s is not far from the conjectured lower limit of 0.08 
for quantum liquids. Moreover, its T-dependence seems 
to indicate a minimum value when approaching Tc from 
above which would resemble rather generic behavior of 
substances in the vicinity of a critical point. This fea- 
ture is also present for the i^-potential, albeit the T- 
dependence is less pronounced; of course, the values are 
also much larger compared to the {/-potential, by about 
a factor of 2-4. In fact, the rj/s value for F at 2Tc is 
rather close to perturbative calculations [i^, Q where 
77/s ~ 1 with little T-dependence, indicating that resum- 
mation effects in the T-matrix do not play a large role 
under these conditions. 



C. Charm-quark number susceptibility 

Quark- number susceptibilities, Xqi which we already 
utilized in connection with the c-quark diffusion coeffi- 
cient, Eq. (P^ . can be computed rather accurately in 
IQCD [4^ and are therefore of great interest to constrain 
effective models of the QGP. For example, HQ number 
susceptibilities have been used to extract effective in- 
medium charm- and bottom-quark mass corrections by 



fitting a zero-width quasiparticle expression, Eq. ([2 
to the IQCD results [43 • Here, we carry out a full off- 
shell calculation including finite-width effects through the 
single c-quark pro pag ators figuring into the charmonium 
spectral function |47| . 



CJO 

'~ tJ 2^ 1 



cM-E/T) 



PooiE,0) . (24) 



For high T ^ rnq, Sq, this quantity reduces to 



XciT) 



2Nr 



6 



(25) 



In Fig. [5] we summarize our results for the HQ sus- 
ceptibility and compare to IQCD data as well as the 
zero- width quasiparticle limit, Eq. (|19p . used in previous 
estimates. When plotted as a function of temperature 
over in-medium c-quark mass (left panel), the results for 
the /7-potential are slightly above the ones for F (note 
that the in-medium c-quark mass is 10-20% larger for the 
t/-potential, mostly due to the larger mass correction re- 
sulting from ?7oo/2 compared to Foo/2). The full results 
using U are significantly above an off-shell calculation 
where the imaginary parts are put to zero, which in turn 
agree very well with the zero-width quasiparticle limit 
(dashed line) . When plotted as a function of temperature 
(right panel in Fig. [5]), it turns out that the results for 
F are somewhat above those for U, due to the smaller 
mc{T). On the other hand, the finite- width effects in 
the full results using U produce an increase in Xc over 
the no-width limit which corresponds to a c-quark mass 
decrease of ca. 150-200 MeV in the zero-width quasipar- 
ticle expression. Thus, finite widths considerably affect 
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the extraction of the in-medium c-quark mass from the 
susceptibihty. This effect also leads to a significant im- 
provement in the comparison to IQCD results: the full 
results with U roughly lie in the uncertainty band encom- 
passed by the IQCD results while the full results with F 
tend to lie at the upper end of that band. At the low- 
est considered temperature, the results for U seem to lie 
below the Nt—8 IQCD computations, but as indicated in 
connection with the large-r limit of the S-wa.ve charmo- 
nium correlator ratios, we believe that further width (and 
coupled-channel) corrections need to be included close to 
Tc before more precise conclusions can be drawn. At least 
the underestimate of the IQCD data in our calculations 
for Xc and the large-r limit of Rps,v is consistent. 



V. CONCLUSIONS 

We have conducted a study of charmonium and open- 
charm properties in the Quark-Gluon Plasma using a 
thermodynamic T-Matrix formalism^. For the first time 
in the context of heavy quarks in the QGP, we have im- 
plemented this scheme selfconsistently at the one- and 
two-body level (i.e., selfenergy and scattering amplitude), 
including microscopically calculated off-shell effects, in 
particular imaginary parts. Within the the Matsubara 
formalism, the two-particle propagator in the T-matrix 
equation automatically generates "zero-mode" contribu- 
tions, i.e., scattering off pre-existing charm quarks in the 
heat bath. Following our earlier work the two-body 
input potential was constructed using a field-theoretical 
model for Coulomb and confining forces with its 4 pa- 
rameters fitted to finite-temperature lattice QCD data 
of the color-averaged heavy-quark free energy. As limit- 
ing cases of the interaction strength we have considered 
the resulting free and internal energy as an underlying 
potential. Once the input potential and the bare charm- 
quark mass are fixed (the latter is adjusted to reproduce 
the charmonium ground-state mass in vacuum) , there are 
no further tunable parameters in the heavy-heavy and 
heavy-light sector of our approach; heavy-quark inter- 
actions with gluons are treated perturbatively but turn 
out to play a minor role (even with 0^=0.4). We are then 
able to comprehensively compute hadronic spectral func- 
tions in different charmonium and D-meson channels, as 
well as heavy-quark selfenergies and transport properties 
in the QGP. Specifically, we have applied these quan- 
tities to calculate euclidean correlator ratios for char- 
monia (consistently including zero-modes in scalar and 
axial/vector channels) and the charm-quark number sus- 
ceptibility, which have both been "measured" with good 
accuracy in thermal lattice QCD. Overall, we find an en- 
couraging degree of agreement of our results with those 
in IQCD, especially at temperatures T > 1.5 Tc. (possi- 



bly with a preference for using U as potential) , keeping in 
mind that neither F nor [/, as computed in IQCD, nec- 
essarily provide a rigorous definition of a two-body po- 
tential, as both quantities are computed from differences 
of thermal averages. Especially close to Tc a more com- 
plete description of charm and charmonia remains a chal- 
lenge, requiring further refinements such as the inclusion 
of (pre-) hadronic states in the QQ and Qq T- matrices 
(coupled channels), a nonperturbative treatment of in- 
teractions with thermal gluons and a careful assessment 
of retardation effects. The large entropy contribution 
in the HQ free energy close to Tc indeed suggests many 
additional states to play a role. Another extension of 
our approach concerns the 3-momentum dependence of 
charmonium correlators/spectral functions, for which in- 
teresting IQCD results are now becoming available [48 1. 

Our analyses and findings reveal the intricate relation 
between charmonium bound-state properties and heavy- 
quark diffusion in the QGP, which, in particular, cannot 
be captured by perturbative treatments. The very same 
(resummed) force that generates resonance-like correla- 
tions (or threshold enhancements) in charmonium spec- 
tral functions - crucial for properly describing the IQCD 
correlator ratios ~ is operative in low-momentum charm- 
quark scattering (closely related to zero modes), and thus 
most likely instrumental in reducing the pertinent ther- 
mal relaxation times as re quir ed in phenomenological ap- 
plications to RHIC data |l4l . |49| . It remains an open 
question how large c-quark momenta need to be for ra- 
diative processes to become competitive with elastic scat- 
tering [iiiii. 

In conclusion, we believe that studying the many-body 
physics of heavy-quark systems in the QGP can indeed 
reveal valuable insights into the medium modifications 
of the basic QCD force via its effects on bound-state 
properties and heavy-flavor transport. The nonpertur- 
bative nature of the problem, especially in the vicinity of 
Tc, requires effective approaches whose predictive power 
greatly benefits from constraints obtained from thermal 
lattice QCD. This will hopefully pave the way toward 
illuminating and quantifying the properties of the QGP 
even in the case of rather strong coupling, and help un- 
derstand some of the fascinating phenomena observed in 
ultrarelativistic heavy-ion collisions. 
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